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Abstract 

We calculate the amplitudes of J PC = 3 meson production in diffractive DIS 
within the /^-factorization approach, with a particular attention paid to the p3(1690) 
meson. We find that at all Q 2 the ,03(1690) production cross section is 2-5 times smaller 
■ than the /}(1700) production cross section, which is assumed to be a pure D-wave state. 

Studying o"l and <jt separately, we observe domination of p3 in ox and domination of 
/3(1700) in gt and offer an explanation of this behavior in simple terms. We also find 
very strong contributions — sometimes even domination — of the s-channel helicity 
violating amplitudes. The typical color dipole sizes probed in p% production are shown 
\ to be larger than those in the ground state p production, and the energy dependence of 

P3 cross section turns out to be much flatter than the p production cross section. All the 
. conclusions about the relative behavior of ,03 (1690) and ,0(1700) mesons are numerically 

Qh| stable against variations of input parameters. 

CD . 

-p. 1 Introduction 

> : 

^ . Diffractive production of vector mesons (VM) in DIS YP ~~ > V"P \Y — P-,4>-,J/' l P e ^ c -) is a 
very active field of research (see recent review d] and references therein). Studying the Q 2 - 
behavior of the VM production cross sections, one can learn about the transition from soft to 
hard regimes in strong interactions, while their energy dependence reveals the Regge properties 
of the Pomeron exchange. The rich set of possible helicity amplitudes 7(A 7 ) — > V(Xv) allows 
one to study the spin properties of the reaction and to test the s-channel helicity conservation 
(SCHC). 

The main focus of this research has been on the production of the grounds state mesons, 
while diffractive production of excited states did not enjoy much attention. Perhaps, the 
most studied case so far was the production of radially excited mesons V(2S). Remarkable 
consequences of the presence of a node in the radial wave function described in [2] were nicely 
confirmed by HI measurements of diffractively produced ip(2S) j^j. 

Similar experimental studies of excited p mesons are expected to be even more rewarding. 
First, diffractive production of excited p mesons probes the dipole cross section at larger dipole 
sizes than the production of ground states. For example, in analysis of [1] dipole sizes up to 
2 fm were important. Such a unique sensitivity of these reactions to soft diffraction can help 
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understand the phenomenon of saturation, which is now a hot topic of debates (see j3] and 
references therein). Another handle offered by diffractive production of p' is their possible 
help in resolving the long standing puzzle of the radial/orbital excitation assignment, as well 
as a possibility of a hybrid component in the p(1450) and p(1700) mesons, see [Hj. 

Diffractive production of excited p' mesons has been observed for a long time in a number of 
fixed target experiments with relatively high energies. Diffractive production of p'(1600) was 
reported in 7r + 7r~ [7] and 4-7T jH] final states (for reanalysis of these data in terms of p(1450) and 
p(1700) mesons and for references to earlier experiments at lower energies, see 0). These states 
were also studied in a recent Fermilab experiment E687 ^H] both in 2ir and 4tt channels. Fi- 
nally, the OMEGA Collaboration has succeeded to measure diffractive photoproduction cross 
section of the p3(1690) (known then as g(1690) meson) via the a2(1320)7r — > t]7t + tt~ diffractive 
final state [TT]. However, all these experiments gave only the value of the photoproduction 
cross section, and no energy dependence, Q 2 dependence, or helicity structure of the reaction 
was studied. This gap was partially closed by the HI measurements of p' electroproduction at 
4 < Q 2 < 50 GeV 2 [12], but due to low statistics the results presented had large errorbars. 

Diffractive production of excited states has not received too much attention also from 
theory. Early theoretical discussions were limited to vector dominance models and its off- 
forward upgraded versions, see j9|. The pQCD based calculations of diffractive production 
of the ground state vector mesons were developed in mid-90 's and were almost immediately 
extended to the case of radial excitations, since the principal effect there is the presence of 
the node in the radial wave function |2j. However, for a long time no microscopic calculation 
of the orbitally excitated vector mesons was available. 

The situation was aggravated by an understanding that production of .D-wave vector 
mesons (as well as J PC = 3 mesons) should be suppressed by Fermi motion, as its radial 
wave function vanishes at the origin. This suppression was believed to be sufficiently strong 
and even prompted the authors of 4J to consider diffractive p(1450) and p(1700) production 
neglecting in both cases the D-wave contributions altogether. Only in ^3] were the .S-wave 
and .D-wave vector meson production amplitudes calculated within the fc t -factorization ap- 
proach, however at that time the absence of convenient parameterizations of the unintegrated 
gluon density — the key input quantity — impeded numerical predictions. 

During the last years, several fits to the unintegrated gluon density appeared jHj, [TK] . 
This allowed for the first estimates jTHJ of the purely D-wave VM production cross sections, 
which showed that at small to moderate Q 2 the production rates of the D-wave and 2S p' 
states are roughly comparable. This was not surprising, since similar conclusion was drawn 
in during the famous splitting the p'(1600) into p(1450) and p(1700). For a more detailed 
analysis of the D-wave vector meson production in fc r factorization, see [T7] . 

The approach developed in jT2] can be applied also to the diffractive production of spin-3 
mesons. The D-wave vector meson and the spin-3 meson can be viewed as spin-orbital splitting 
partners and can be described within the same formalism. The only modification required to 
proceed from D-wave VM to spin-3 meson is that of the qq coupling to the final meson: 

uv» u -v;^uC^u-t; vp , (i) 

where and T^ vp are the polarization vector for spin-1 and polarization tensor for spin-3 
mesons, respectively. The spinorial structure T>^ determined in corresponds to the pure 
D-wave vector meson, while the structure C fJ,up to be found should correspond to D-wave 
spin-3 meson. 
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In this paper we report the first microscopic derivation of the J PC = 3 production 
amplitudes within the /^-factorization approach. We focus on the p3(1690) production and give 
predictions on the Q 2 and W dependence of the cross sections, on the <7l/<7t decomposition, 
and on the role of s-channel helicity violating amplitudes. We also compare the p3(1690) 
cross sections with those of p(1700), which is assumed to be purely D-wave vector meson, and 
observe a number of remarkable distinctions. These should prove useful in disentangling these 
two mesons in experiment, especially in the case of low statistics. 

The paper is organized as follows. In Section 2 we show how spin-3 mesons are described 
and present analytic expressions of the spin-3 meson production amplitudes. In Section 3 we 
perform twist analysis of the amplitudes in the forward case, compare the results with those 
of D-wave vector mesons, and discuss the effects of s-channel helicity violating amplitudes, 
which will appear in the non-forward scattering. In Section 4 we present numerical results for 
P3 and compare them with the corresponding ground state p and p(1700) cross sections. In 
Section 5 we discuss the typical color dipole sizes probed in p 3 production, compare our results 
with experimental data available and comment on future possibilities. Finally, in Section 6 
we draw our conclusions. 

2 Amplitudes of spin-3 meson production 
2.1 Kinematics and notation 

We use the usual notation for kinematical variables. Q 2 is the photon's virtuality, W is the 
total center-of-mass energy of the 7*p collision. The momentum transfer from proton to photon 
is denoted by A M and at high energies is almost purely transverse: —A 2 = \t\ ~ \t'\ = A 2 . The 
transverse vectors (orthogonal to the 7*p collision axis) will be always labelled by an arrow. 

Diffractive production of meson V with mass my can be treated in the lowest Fock state 
approximation as production of the corresponding qq pair of invariant mass M ^ my, which is 
then projected, at the amplitude level, onto the final state. Within the leading log- accuracy 
the higher Fock states are reabsorbed into the evolution of the unintegrated gluon density 
(or color dipole cross section). A typical diagram to be calculated (see [I]) contains the 
valence quark loop, with integration over the quark transverse momentum k and its fraction 
of photon's lightcone momentum z, and the uppermost gluon loop, with the integration over 
transverse momentum k. A convenient choice is to assign momentum k + zA to the quark and 
—k + (1 — z)A to the antiquark, which ensures that even at non-zero A the qq invariant mass 
is M 2 = ^|+^ _ it j s a i so convenient to consider the relative qq momentum p^ = {k q — kq)^/2 

in the qq pair rest frame, where it reduces to the 3-momentum p = {k, k z ) with k z = ^—M. 
Such 3-dimensional vectors, which always refer to the qq rest frame, will be given in bold. 

Finally, throughout the text the ground state vector mesons (which will be always under- 
stood as IS* states) will be generically labelled by V or Vis, the pure D-wave vector mesons 
will be labelled by Vd, while the J PC = 3 states of the same quarkonium will be marked 
as V3. When we speak of D-wave mesons, we will always assume D-wave vector mesons, not 
spin-3 mesons (although in V3 the qq pair also sits in the D-wave). In order to avoid excessive 
subscripts, my will refer to the mass of the meson being discussed. In particular, for the p 
system we will speak of p 3 and po- The former refers to the physical ps(1690) state, while 
the identification of the latter with p(1700) is clearly model dependent and is done only for 
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purposes of comparison. The effects of S'-wave admixture in p(1700) are considered in |17J. 
2.2 Description of a spin-3 meson 

A spin-3 particle is described with the rank-3 polarization tensor T Ml ' p , which must be symmet- 
ric and traceless in any pair of Lorentz indices. States with any given helicity can be written 
as simple combinations of well-known polarization vectors e^, A = +1, 0, —1 (see below more 
on eft): 

rrlivp _ P p V P rpfivp _ f U v p\ 

m p = ^= (2{e» + ey } + {e^ + e p _}) , = ^= (2e^eg + {e^e* }) , (2) 

where curly brackets denote symmetrization. 

To construct the spinorial structure C pI/p , recall first that both quark and antiquark spinors 
can be treated on-mass-shell (see details in ^El)- Therefore, C^ up can be contructed from two 
independent structures: ^p v pp and pPp v pP (symmetrization over indices is assumed). When 
constructing C tiup , we should make sure that it represents purely L = 2 state without admixture 
of L = 4 wave, which is in complete analogy with construction of S and .D-wave vector mesons 
performed in [T^J. Note that we explicitly rely on the lowest Fock state approximation of the 
meson; if higher Fock states are taken into account, such simple picture is lost. 

qwp j n (JJJ can k e constructed in a most trasnparent way in the non-relativistic case. 
Instead of the D-wave vector meson structure (p' D^a^ (p, where D l i = ?>p l pp — 5 lJ p 2 , we have 
now </r D % io k tp. Recall that this structure will be contracted with polarization tensor T^ k so 
that one should not worry about symmetrization. Since this polarization tensor is traceless, 
the coupling simplifies to ip^ a k ip p l p> , which is just a tensor product of the S- wave coupling 
for vector meson and the p % pi term. We do not have to keep track of the overall normalization 
factors, since they can be always absorbed in the definition of the radial wave function. 

Returning to the fully relativistic case, we can now write the spinorial structure for spin-3 
meson coupling 

u u ■ = u <S W u-T^ p = uS^u-r^ r^ = T FpP y . (3) 

As the last form of this coupling shows, the spin-3 amplitudes are easy to construct once we 
know the amplitudes for the S'-wave vector meson. 

Given the spinorial structure, one can now calculate various amplitudes with spin-3 meson. 
In such calculations, the radial wave function ip(p 2 ) will appear, whose normalization is 

Apart from coefficient 1/15, this normalization condition coincides with the .D-wave vector 
meson normalization condition. 

We underline that our approach to describing polarization states of the final meson is 
explicitly rotationally- invariant. We use identical radial wave functions for all the polarization 
states of the final meson, and in this way the normalization condition (J3J holds for an arbitrary 
polarization state of the spin-3 meson. An important part of the rotation-invariant description 
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is that the transversity condition must be imposed at the level of the qq pair. This leads to the 
concept of the running longitudinal polarization vector eg(M), such that it is orthogonal to the 
4-momentum of the on-mass shell qq pair with invariant mass M rather than the momentum 
of the final meson. The calculations with the fixed longitudinal polarization vector often found 
in literature break the rotation invariance. In technical terms, the fixed polarization vector 
leads to a mixing of the longitudinal spin-1 state and spin-0 states. 

Potential models suggest that the radial wave functions of the spin-orbital partners should 
be very similar. One can assume, as a starting approximation, that their shapes are identical. 
This assumption leads to 

^ 3 (p 2 ) = v / T5^d(p 2 ), (5) 
which will be useful for comparison of spin-3 and .D-wave vector meson production. 



2.3 Generic amplitudes 

A generic form of the helicity amplitudes 7*(A 7 
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Here cy is the flavor- dependent average charge of the quark, the argument of the strong cou- 
pling constant a s is max[z(l — z)(Q 2 + M 2 ), k, 2 }, and J-(xi, x 2 , k, A) is the skewed unintegrated 
gluon distribution, with x\ ^ x 2 being the fractions of the proton's momentum carried by the 
uppermost gluons. The appearance of skewed (or generalized) parton distributions is charac- 
teristic for scattering processes that change the mass/virtuality of the projectile JTHj , see also 
recent reviews [TH]. In the case of meson production, its use is important due to x 2 <C X\ 
and has been incorporated in the collinear factorization approach [T%1 |2*U] as well as in the 
factorization approach with non-zero transverse momenta of quarks taken into account fZT\ . 
In the /^-factorization approach, the skewness is transferred to the unintegrated distributions. 

The integrands for spin-3 mesons I\ s .\ can be written in terms of the corresponding 
integrands for vector mesons: 
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Note that both p M and e x depend on A, still this dependence of their scalar product vanishes 
due to the Lorentz invariance. The integrands I^-x are: 
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The integrands for helicity —1 can be obtained from those with +1 by replacement of k + — ^ k_ 
and $ 1+ — > $!_ (with no extra minus sign that would appear only at the level of amplitudes!). 
Here function $ 2 describes transition of virtual photon into the qq states with \ q + Xq = A 7 * , 
whereas $i describes transition of transverse photons into the qq states with X q + Xq = 0, in 
which the helicity of the photon is carried by the orbital angular momentum in the qq state: 
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where r = k + (2z - l)A/2 and Q = z(l - z)Q 2 + m 2 



3 Large Q 2 , my analysis 

The above expressions can be integrated numerically. However, before describing these results 
it is useful to study analytically the case where both Q 2 and my are much larger than any 
soft scale, while Q 2 /my can be arbitrary. In this approximation one expands the hard scale 
Q 2 around Qq = \{Q 2 + fn v ), as well as performs expansion in powers of small Fermi motion 
of the qq pair, k 2 zl k 2 -C my. We will call this approximation the "twist" expansion. We start 
with the forward case, A = 0, where only the s-channel helicity conserving amplitudes with 
A 3 = A 7 survive, and find the ratio ctl/vt as well as relation between V 3 and Vp production 
cross sections. After this, we discuss the role of helicity violating amplitudes. 

3.1 Twist expansion for the forward case 

We consider the two non-zero integrands in this approximation Q and (jlUJI and note that after 
dflp angular averaging all terms in each of these expressions give comparable contributions, 
cx P 4 k 2 /Qq, differing only by numerical coefficients: 
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It is curious to note that the leading-twist contribution in the transverse amplitudes vanishes, 
and one is left with the subleading term. This cancellation does not occur in the longitudinal 
amplitude, which leads to un abnormally large value of the ratio (Jl/o't- 

b -<"> <- 27 (i+ Q2 V»i an 

Rlt = ^-Q^-Y{ 1+ 3K) >>l - 

This must be confronted with R^j? = 1 for the ground state mesons and, even more remarkably, 
with Rlt <C 1 for D-wave vector mesons, evaluated within the same approximation. 

We stress that such a peculiar Q 2 -dependence of the ratio ul/ct (H3) is entirely due to the 
heavy-meson approximation we used. Allowing for the longitudinal quark motion will restore 
the leading twist contribution to the transverse amplitude. What is expected to remain, 
however, is the overall smallness of the transverse amplitude and, therefore, a large numerical 
value of Rlt- 



3.2 Spin-3 vs. D-wacve vector meson 

Within the twist expansion, the Vp and V 3 production amplitudes are proportional to J d 3 p p"Vz)(p 2 ) 
and / G? 3 pp 4 ^3(p 2 ), respectively. Assuming (JHJ), one can relate the Vd and V3 production cross 
sections. The results for the longitudinal and transverse cross sections, separately, are 
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where we assumed the masses of the two states to be equal. One sees that the longitudinal 
cross section is dominated by V3 meson, while the transverse one is dominated by Vd. In some 
sense, these two mesons "mirror" each other: where V3 is suppressed, Vd dominates and vice 
versa. 

This "mirror" behavior can be in fact understood in simple terms. Consider, for instance, 
the T — > T transition. Note first that although the integrands for Vd derived in flB] look 
very differently from those of V3, they can be written in a form similar to ([TJ: uD^u ■ V^ = 



uS^u ■ D^p ■ V p , where ■ Vl 
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The key point is the opposite signs in front of the second term. It turns out that the contribu- 
tions of all three terms in are of the same sign and of the same order of magnitude, so that 
they interfere constructively in Vd production. In the case of V3, they interfere destructively, 
which leads to suppressed 

For the longitudinal amplitude, the similar change of signs strongly suppresses the result for 
D-wave mesons, enhancing it in the spin-3 case. We see that there are good reasons to expect 
such a "mirror" behavior of V3 and Vd just on the basis of their spin-angular composition. 
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3.3 The role of s-channel helicity violation 

The approximate conservation of the s-channel helicity in diffractive reactions is due to two 
reasons. First, in diffraction the helicity properties of the target and projectile are uncorre- 
lated, and in the forward case strict SCHC holds separately for the projectile and the target. 
At A 7^ the s-channel helicity non- conserving (SCHNC) amplitudes depend on the momen- 
tum transfer as |A|' A7 ~ Al/ ', which makes them small within diffractive cone. Second, at high 
energy the helicity is conserved at the parton level, so in order to produce helicity flip the 
transverse motion of constituent must come into play. This produces extra factors like k 2 /M 2 , 
further suppressing helicity violation, especially for heavy quarks. Nevertheless, small viola- 
tion of SCHC has been observed at HERA in the case of light vector mesons; at the amplitude 
level, its relative magnitude was estimated to be ~ 10% |22] . 

In the case of V3 production, the effect of SCHNC must be more important, just as it was 
for D-wave vector meson production The integration of the quadrupole term kills the 
leading contribution to the SCHC amplitudes (fT^j) and (|T3j). and the SCHC amplitudes get the 
same suppression due to Fermi motion as SCHNC ones. Moreover, partial cancellation among 
several terms discussed above suppresses the T — * T amplitude, while the helicity violating 
amplitudes do not suffer such cancellation. Finally, pay attention to numerical factors like 
l/\/l5 in the amplitude A +1 . +1 , which are absent, for example, in the A +3 . + i amplitude and 
take into account the large number of various helicity violating amplitudes for spin-3 meson 
production. 

Thus, one can anticipate that the helicity violating amplitudes can generate a significant 
portion of the overall cross section. One should not even be surprised to see them dominate 
in the transverse cross section, especially at small Q 2 . Therefore, the above twist analysis is 
meant only to guide the eye and should not be used for quantitative discussion. 

4 Numerical study 

In this section we present numerical results for the particular case of p% production. 
4.1 Input 

In order to integrate (jBJ) numerically, one needs to specify models for the unintegrated gluon 
density and the meson wave function. We related the skewed unintegrated gluon density with 
non-zero momentum transfer to the forward unintegrated gluon density by 

t x a , k + 1a, -a + \ a) = t (o.4i^±^, a) exp , (20) 

where 631P includes contributions from two-gluon formfactor of the proton and from the effec- 
tive Pomeron trajectory, as described in detail in 1J. Although fixing the exact numerical value 
of the shift coefficient (0.41) is beyond the log^ accuracy, its introduction is phenomenologically 
motivated and, after all, it can be viewed as yet another parameter in our parametrization. We 
did not try varying this parameter to obtain a better fit. The parametrizations for the forward 
unintegrated gluon density were borrowed from [TSJ. Note that the /^-factorization approach 
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itself does not require Q 2 to be large, and in the soft region, Q 2 < 1 GeV 2 , the words "unin- 
tegrated gluon density" should be understood simply as an appropriately normalized Fourier 
transform of the dipole cross section. 

As for the radial wave function, we choose the simple Gaussian Ansatz WF normalized 
according to (J4*J) . with the only free parameter, the "size" a of the wave function. In contrast to 
the vector meson case, dileptonic decay p 3 — > e + e~ cannot proceed via one-photon annihilation 
and does not help us fix a. However, appealing to the argument that the radial wave functions 
of spin-orbital partners should be similar, we can take a 3 = od, the latter being extracted 
from r(p(1700) -> e+e"). 

4.2 Level of accuracy anticipated 

Our experience with diffractive production of ground state vector mesons within the same 
approach tells us that variation of the input parameters changes the absolute values of the 
cross sections by a factor of < 1.5, while the accuracy for the observables that depend on the 
ratios of the amplitudes is even better PES]. The principal source of uncertainty was found 
to be the final meson wave function, especially its density near the origin. The sensitivity of 
the results to particular parametrizations of the unintegrated gluon density presented in jTH] 
(from both DGD2000 and DGD2002 sets of parametrizations) was found to be weak. 

In the present case, the results are expected to be less stable with variation of input due 
to the presence of various cancellations. The main source of instability is the poorly known 
value of the pr> dileptonic decay width. The data available give T(p D — > e + e~) ~ 0.1-0.6 
keV (assuming that p(1700) is indeed D-wave qq state). The possibility that p(1700) has 
significant contributions from radially excited qq and from possible hydrid state, as well as 
taking into account extremely large NLO corrections j2Hj for this decay, make the situation 
even less definite. 

All the curves to be presented below were calculated for T(pr> — > e + e~) = 0.14 keV. This 
value corresponds to the value of r(e + e~) • Br{jt + -n~ -) = 29±i| eV obtained in [21]. In order to 
see the effect of this input parameter, we calculated cross sections both for p 3 and po for the 
dileptonic decay width in the interval 0.14 — 0.7 keV. Increasing T(p£, — > e + e~), we observed 
some suppression of the cross sections at small Q 2 and their significant growth at Q 2 > 1 
GeV 2 , especially in the case of a-r- The effect is strong, and we conclude that the numerical 
results for the absolute values of the cross sections are trustable only within factors of ~ 2-3. 

We stress, however, that variation of the input parameters produced absolutely the same 
shifts in the po production cross sections. This should be expected, because the relation 
between the two mesons is dictated primarily by the similarity of their radial distributions 
and by spin-angular relations of type of (|19j) . These relations are essentially insensitive to 
details of the model, as long as we treat p(1700) as a predominantly qq pair in the .D-wave 
state. We conclude therefore that the numerical values of our predictions for ratios between 
P3 and pe> are more stable, approximately within a factor of 1.5-2. 

4.3 Q 2 and ^-dependence 

We calculated all the helicity amplitudes (jTjl- lfTUI) for spin-3 meson and compared its production 
rate with that of pr, and pis- All the cross sections are calculated at W = 75 GeV and 
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are obtained from numerical integration of the differential cross sections within the region 
< \t\ < 1.05 GeV 2 . 

Fig. [T] shows the ratios of the excited to ground state meson cross sections cr(p 3 )/a(pis) 
and cr(pu) j 'c(pxs)- Both ratios are an order of magnitude smaller than unity, and the po cross 
section is noticeably larger than that of p%, especially at Q 2 ~ 1 GeV 2 . At larger Q 2 , the ratio 
a(pD) I &(pz) ~ 2. Thus, if one intends to extract p(1700) properties from diffractively produced 
multipion states around invariant mass ~ 1700 MeV, one cannot neglect contamination by 
the p 3 state. 
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Figure 1: Predictions for the ratios of p 3 (solid line) and po (dashed line) to p production 
cross sections as a function of Q 2 . 

Difference between pr> and p 3 is better seen if one studies separately longitudinal and 
transverse cross sections, shown in Fig. |21 Here, the solid and dash-dotted lines represent 
the p 3 and p D cross sections, respectively. In the case of p 3 we showed also with the dashed 
lines the contributions of the SCHC amplitudes only. One clearly sees the domination of 
helicity violating amplitudes at small Q 2 in p 3 production. One can even state, on the basis 
of our calculations, that p 3 production at small Q 2 probes diffraction in the regime of strong 
s-channel helicity violation. In the case of longitudinal cross section, the contribution of 
SCHNC transitions becomes small at Q 2 > 1 GeV 2 , since all such amplitudes are of higher 
twist. Helicity violation remains strong for transverse photons even at large Q 2 . 

As mentioned above, in the case of D-wave vector mesons one expects suppression of <7l 
but not in ot- Indeed, our calculations show the domination of the p 3 over pn in o~l at small 
Q 2 <1 GeV 2 , while in ot the po cross section is noticeably larger than p 3 everywhere. This 
is in a qualitative agreement with twist analysis result (J18)) . 

Such a different behavior of prj and p 3 can be seen also in plots of ratio a^l^T as a function 
of Q 2 , shown in Fig. El Here, on the left plot, we showed this ratio for p 3 and po, while on the 
right plot, we showed reduced ratios Rlt = — • jfe, where rut is the mass of the corresponding 
meson. Rlt is small for po and relatively large for p 3 , as was expected from the twist analysis 
(JT7J). We note that the region Q 2 ~ 1 GeV 2 is particularly suitable for distinguishing among 
various p states. 
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Figure 2: Predictions for the longitudinal (left plot) and transverse (right plot) cross sections 
of p 3 (solid lines) and po (dash-dotted lines) production. The contribution to p 3 from SCHC 
amplitudes only is shown with dashed lines. 




Q 2 (GeV 2 ) 

Figure 3: Predictions for the ratios of the longitudinal to transverse cross sections of p 3 (solid 
line) and po (dashed line) production. The left plot shows R = gl/o~t, while the right plot 
shows reduced ratio Rlt = olI^t ■ fny/Q 2 . 



The role of helicity violating amplitudes can be seen also in the | ^-distributions shown in 
Fig. |U for Q 2 = 1 GeV 2 separately for longitudinal and transverse cross sections. The dashed 
and dotted lines show the SCHC contributions only, while the solid lines show their sums. As 
can be expected, the large contributions of SCHNC amplitudes come from the entire i-interval 
shown, while the SCHC amplitudes are strong only within a narrow forward cone. 
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Figure 4: Differential cross sections do^jd\t\ (left plot) and d<j T jd\t\ (right plot) for p 3 at 
Q 2 = 1 GeV 2 as functions of \t\. The dashed and dotted lines show contributions from helicity 
conserving and helicity violating amplitudes, respectively; the solid lines show their sum. 



4.4 Energy dependence 

We checked also the energy dependence of the p% production cross section, which we parametrized 
with a simple power law cr(p 3 ) oc W s . General expectations, driven by the ground state pro- 
duction experience pQ, are that at small Q 2 this exponent should be a small number, and one 
speaks usually of "soft Pomeron" , while in the presence of hard scale it should grow up to 
8 ~ 1, corresponding to the "hard Pomeron". 

For a particular example, we used HERA kinematics and studied the energy behavior of 
the cross section within the range W = 50-200 GeV. We found that at small-to-moderate 
Q 2 , the P3 production cross section (both longitudinal and transverse) slightly decreases with 
energy rise, with typical 5 ~ —0.1 to —0.2. It is only at Q 2 > 5 GeV 2 that 5 becomes positive, 
and it is always smaller than the value of 5 for the corresponding ground state vector meson. 

This decrease is naturally understood in the Regge picture of the Pomeron exchange. The 
differential cross section at non-zero t behaves roughly as 

^ oc W 5 ^ ■ 5(t) = 4[ ajP (t) - 1] « 4[o„(0) - 1 - a' eff ■ \t\)] . 

The value of the effective Pomeron intercept ajp(0) depends on Q 2 and comes directly from the 
parametrizations of the forward unintegrated gluon density, see [15] . At small Q 2 it is about 
«ip(0) — 1 ~ 0.08, and starts noticeably growing only at Q 2 > 2-3 GeV 2 . The effective slope 
of the Pomeron trajectory is ol e ^ w 0.12 GeV -2 with very marginal Q 2 dependence. (Note 
that it differs from the fixed input parameter a' e ^ ~ 0.25 GeV -2 used in our calculations due 
to ant i- shrinkage effects discussed in detail in [25].) Thus, the effective exponent of the energy 
dependence of the integrated cross section 

tfw4[a I .(0)-l-a , e// .<|t|>] 
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is governed not only by aip(O), but also by the typical momentum transfers (|t|) involved. 

In the production of ground state vector mesons dominant contribution comes from SCHC 
helicity amplitudes, which are concentrated within the forward cone \t\ < 0.1 GeV 2 . In the 
present Fig. |U shows vividly, the range of important values of \t\ spans up to 0.5-1 

GeV 2 . One sees that due to such high values of \t\ involved the energy rise exponent S at small 
Q 2 can easily become negative. 



5 Discussion 

5.1 What dipole sizes are probed in photoproduction? 

The 5-wave qq state is naturally orthogonal to the D-wave qq state. For example, if one 
attempts to calculate the inelastic ps — > Pd formfactor at zero momentum transfer, one finds 

m oc / -^^d 2 kr D (p 2 )Mp 2 mi-k 2 ) = J ^d 3 P r D (p 2 )Mp 2 )(2k 2 z -k 2 ) = o. 

The presence of the quadrupole combination 2k 2 — k 2 makes the amplitude zero, as long as 
all other factors under the integral are spherically symmetric. 

In the photoproduction, the wave function of the initial photon is not spherically sym- 
metric, but one can still appeal to the vector dominance model arguments and rewrite the p% 
photoproduction amplitude as 

(7|a(r)|p 3 )oc^ p .(p| ( 7(r)|p3>. (21) 

One can suspect that a similar orthogonality should be at work here, when one considers a 
forward SCHC amplitude at large qq dipole sizes, where the dipole cross section a(r) — > const 
so that all other factors seemingly become spherically symmetric. If this were the case, it 
would mean that the p% photoproduction receives little contribution from large dipoles and 
is a "harder" process than the p photoproduction. However, this is not the case. The most 
essential difference between the pp — > p$p amplitude and (jUJ) is the replacement of the photon 
wave function: 



$0 -> *r 



zl 



2Mk 2 z + k 2 ) - Mkl + (k + k) 2 ) - Mk 2 z + (k - nf) , (22) 



together with a similar replacement of $i. For clarity, we explicitly presented the spherically 
symmetric quantity p 2 as k 2 + k 2 . At small values of k 2 , which correspond to large dipole 
sizes, one obtains: 

M/ 2 ^-2/? 2 [^(p 2 )+^(p 2 )] , (23) 

where derivatives of the radial wave function are taken in respect to p 2 . 

The result (|23|) explicitly lacks spherical symmetry. Thus, even if one uses the leading twist 
contributions, pretending that the p — > p 3 transition is well approximated by non-relativistic 
expressions, one still gets 

/ d 3 p [^(p 2 ) + k 2 ^(p 2 )} ^ ■ (2k 2 z -k 2 )^0. (24) 
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This result is rather natural. The characteristic feature of the high-energy collision is presence 
of a preferred direction: that of the proton's momentum in the vector meson rest frame. 
The transverse and longitudinal dynamics of the quark loop now differ, and this leads, in 
particular, to circulation of purely transverse momentum k in the loop, see (J23|) . which breaks 
the spherical symmetry. Taking into account the Fermi motion makes the expression to be 
integrated even less symmetric. 



W t (k 2 )/W t (0) 




k 2 (GeV 2 ) 

Figure 5: The mapping functions Wt{k 2 ) normalized to the k 2 = point for the p% (solid 
lines) and pis (dotted line) photoproduction. The shaded area shows the distribution of the 
results due to the decay width variation T(pr> — > e + e~) = 0.14-0.7 keV. 

It is useful to check explicitly that large dipoles indeed contribute substantially to p^ 
production amplitude. In order to test this numerically, we change the order of integration in 
(jHI) and represent the forward T — > T amplitude as 

1 r dx 2 

—Im A +1 . +1 = J -pH^U X2, «, 0) • W t (k 2 ) . (25) 

One expects Wt(k 2 ) to have the "smoothed step function" shape: it should be approximately 
constant at small k 2 up to some value Qt , and should decrease quickly as k 2 passes this value, 
see for details of this analysis for the ground state p production. This function "cuts out" 
the important range of gluon momenta, and determines thus the important range of the color 
dipole sizes. The effect of orthogonality — if it were present — would appear in Wt( 

Ki 1 clS ct 

small-/? 2 suppression. 

In Fig.Elwe show Wt^ 2 ) normalized to the k 2 = point for the p 3 and pis photoproduction 
with solid and dotted lines, respectively. The shaded area corresponds to scattering of the 
results due to variation of the decay width T(p D — > e + e~) = 0.14-0.7 keV. We observe no 
small-/? 2 suppression for p%. The k 2 value where W(k 2 )/W(0) hits 1/2 is k 2 ~ 0.1-0.2 GeV 2 
for the p3, and is noticeably smaller than the corresponding value of 0.27 GeV 2 in the pis 
production. One sees that the p% production is indeed a softer process, and the typical dipole 
sizes probed are ~ 1.2 4- 1.5 time larger than in the ground state VM production. 
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5.2 Comparison with experimental data available 

The OMEGA Collaboration at CERN measured the cross section of diffractive photoproduc- 
tion of ps(1690) (known then as g(1690) meson) via the a2(1320)7r subsample of the rpT + 7r~ 
diffractive final state events [TT]. The cross section of 7p — > p3(1690)p — > a 2 (1230)7rp was 
found to be 97 ± 28 ± 21 nb, which allows one to roughly estimate the p3 production cross 
section as a{pfp — > ps(1690)p) ~ 200-300 nb. This result is about 5-10 times below our pho- 
toproduction predictions, which we think is not very bad discrepancy, taking into account 
expected level of accuracy in the soft region. Indeed, most of this cross section we predicte to 
be due to SCHNC, especially the double-flip, transitions. Its magnitude in the soft region was 
predicted by our calculation to be rather large even for ground state vector mesons, but so far 
has been poorly known from experiment. We think that upon understanding better the role 
of SCHNC at small Q 2 with the aid of modern experiments, we can improved the accuracy of 
our predictions. We expect, however, that our conclusion of the strong violation of s-channel 
helicity violation probed in p 3 will survive such an upgrade. In addition, our cross sections 
correspond to integration within < \t\ < 1.05 GeV 2 ; the results will change noticeably if one 
selects another t-interval. 

The OMEGA Collaboration also measured the photoproduction cross section of p' . The 
original data were reanalysed in terms of p(1450) and p(1700) separately in c 9j yielding 
cr(p(1700)) ~ 500 nb. Thus, our result cr(p(1700))/cr(p3) ~ 3 at the photoproduction limit is 
roughly consistent with experiment. Finally, comparing the p3 and the ground state p pho- 
toproduction cross sections, we note that our result a^/cris ~ 0.1 is again not very far from 
experimental value of 0.02-0.03. 

5.3 Comments on experimental possibilities 

The experimental analysis of diffractive production of spin-3 resonances, and in particular, the 
strategy of p3(1690)/p(1700) separation, will depend on the statistics available and the final 
state chosen. 

Should one have the luxury of high statistics, one can do the partial wave analysis or select 
some particular final states, in which one of the two states would dominate. An example of this 
approach is just the OMEGA Collaboration observation of the P3 in the ci2(1320)7r — > i]7i + ii^ 
final state. If the statistics does not allow for such angular dependence or final state analysis, 
one then should look for distinctions in the production of these mesons. In the view of our 
results, it is tempting to make use of ratios (Tl/ct, which are dramatically different for p3 and 
Pd, especially in the small-to- moderate Q 2 region. 

One possibility to separate and ot is given by the Rosenbluth method. It will require 
several runs at different lepton beam energies and might seem impractical at high energies. The 
second possibility could be to do a baby-version of PWA and to study angular correlations 
in final state hadrons. For example, if both mesons discussed are observed in 7r + 7r~ (the 
corresponding branching ratios are not dominant, but still sizable), one could study the single- 
differential angular distribution W[co^ff). This measurement will give spin density matrix 
element Tqq, from which one recovers a^joT- Alternatively, one can search for a similarly 
revealing angular dependence in 4ir final states, the dominant decay channel of both mesons. 

Another issue, which requires taking into account the p3 meson, is the recent observation 
of a narrow dip structure in diffractively photoproduced 37r + 37r~ states at M 67r « 1.9 GeV in 
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Fermilab E687 experiment |2Zj- Although the detailed mechanism of its appearance remains 

unsettled, the very recent analysis [2E| sees it as a result of interplay of several resonances 
with J PC = 1 (including p(1700)) and a background. This analysis was explicitly based 
on the vector dominance idea and explicitly uses the assumption that the 67r spectra in e + e~ 
annihilation and in diffractive photoproduction are essentially the same (apart from kinemat- 
ical factors). The results presented here clearly show that this is a risky assumption. The p 3 
meson does not couple to the single virtual photon, yet it should be produced diffractively at 
a rate comparable to that of p(1700). Although it cannot produce any interference pattern 
with J = 1 states, its own contribution can affect the results of a very delicate analysis of [28 . 




6 Conclusions 

We calculated the cross section of the exclusive production of J PC = 3 mesons in diffractive 
DIS within the /c t -factorization approach. The results were compared with the cross section 
of the D-wave state vector meson of the same quarkonium. We exemplified the general ex- 
pressions with a detailed numerical study of the p system, where the p 3 (1690) state is almost 
degenerate with the p(1700) meson, whose structure is arguably dominated by the qq pair in 
the D wave. 

The absolute values of the cross sections suffer from uncertainties of the input parameters, 
in particular, of the pn — > e + e~ decay width, and we can be sure only in the order of magnitude 
of these results. However, in what concerns the relative production rates of p 3 and pz>, our 
conclusions are much more certain. Our results allow us to formulate the following predictions, 
which are stable against variations of the model parameters: 

• In typical HERA kinematics, the ratio of production cross sections taken at equal Q 2 is 
ct(pd) I 'v(pz) ~ 3-5 at small Q 2 , decreasing to ss 2 at larger Q 2 . Thus, when extracting 
the properties of p(1700) from multipion final states, one cannot simply neglect the p 3 
contribution. 

• p 3 and pd show completely different patterns in cjl-vt decomposition: p 3 dominates in 
the longitudinal cross section, while pr, dominates in the transverse cross section. The 
ratios R = cr L /(jT for p 3 and pr> differ by more than one order of magnitude. This 
dramatic difference can be traced back to spin-angular properties of these two mesons, 
see (USD - 

• The role of s-channel helicity violating amplitudes is extremely important, especially 
in the transverse cross section. At small-to-moderate Q 2 the helicity violating ampli- 
tudes even dominate over the SCHC ones. Thus, production of p 3 offers an interesting 
possibility to study diffraction in the regime of strong s-channel helicity violation. 

• Due to large color dipole sizes probed in the p 3 production and large region of relevant 
momentum transfers, \t\ < 1 GeV 2 , the energy dependence of p 3 production cross section 
is less steep than in the case of p. At small Q 2 , one might even observe decrease of p 3 
cross section with energy rise. 

We find no surprise in numerical stability of the first two conclusions, since they are essentially 
driven by very basic relations: similarity of radial wave functions for p 3 and po, the spin- 
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angular composition of these mesons, see (|19|). and quadrupole suppression of the leading 
contributions in SCHC amplitudes (fT2"j) and (fT3*j) . 

In addition, confronting our predictions for photoproduction with the fixed target data 
available and observing them to agree within the anticipated accuracy inspires hope that we 
grasp the essential physics of this reaction in our approach. We are looking forward to seeing 
experimental checks of our predictions. 

It is our pleasure to acknowledge fruitful discussions with Etienne Burtin, Roberto Fiore, 
Nicole d'Hose, Dima Ivanov, Kolya Nikolaev, Alessandro Papa, Andrzej Sandacz. The work 
of I.P.I, is supported by the INFN Fellowship, and partly by INTAS and grants RFBR 05-02- 
16211 and NSh-2339.2003.2. 
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